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SUMMARY 


Higher order spinning modes must be considered in the design of efficient 
noise suppressors with outer wall treatment such as in an engine inlet. These 
modes are difficult to measure and are in fact Impossible to resolve with flush 
mounted wall microphones. An alternative liner design procedure is presented 
here which potentially circumvents the problems of resolution in modal measure- 
ment. The method is based on the fact that the modal optimum impedance and the 
maximum possible sound power attenuation at this optimum can be expressed as 
functions of cutoff ratio alone. Modes with similar cutoff ratios propagate 
similarly in the duct and in addition propagate similarly to the far field. 

Thus there is no need to determine the acoustic power carried by these modes in- 
dividually, and they can be grouped together as one entity. With the optlmimi 
impedance and maximum attenuation specified as functions of cutoff ratio, the 
off-optimum liner performance can be estimated using a previously published 
approximate attenuation equation. 


INTRODUCTION 


The need to consider higher order spinning modes in the design of aircraft 
inlet suppressors with wall treatment only has been demonstrated in references 
1 and 2. Using spinning modes in the propagation theory to simulate an engine 
inlet requires information on the modal power distribution, which is very diffi- 
cult to measure. Assumptions of equal modal amplitude (ref, 3) or equal modal 
power (refs. 1 and 4) have been made. These assumptions may be valid for static 
test data (ref. 5) where the dominant source of noise may be from the inter- 
action of the rotor with random inflow disturbances. However, in flight the 
character of the noise source changes considerably (ref. 6) and the modal struc- 
ture giving valid liner designs has yet to be established. 

Because of the difficulty of modal measurement, an alternative and more 
easily used method has been proposed (ref. 7). This method involves the use of 
the distribution of acoustic power as a function of mode cutoff ratio (hereafter 
called acoustic power-5 distribution) rather than the actual modal power distri- 
bution itself. This is much simpler since many modes may have nearly the same 
cutoff ratio and need not be separated because they all behave the same when 
liner design is considered. This similar behavior is demonstrated by showing 
that the optimxim wall impedance and the maximum possible sound power attenuation 
obtained at this optimum can be expressed as functions of cutoff ratio alone. 
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When these quantities are expressed in this way there is only a very small re- 
sidual modal dependence (lobe m and radial mode y), which can be ignored in 
a multimodal liner design. 

It was established in referenee 7 that modal optimum wall impedance was 
intimately related to cutoff ratio. The reference also implied that maximum 
attenuation and the radiation pattern were dependent upon cutoff ratio. In this 
paper the method will be developed into a quantitative tool useable for liner 
design. Approximate expressions are provided in terms of cutoff ratio for the 
optimum impedance and the maximum possible attenuation. Since no explicit modal 
identity occurs in these required inputs, modal decomposition of the noise 
source is replaced by the acoustic power-^ distribution which treats all simi- 
larly propagating modes as a single entity. These equations also contain the 
usual input quantities such as flow Mach ntimber, boundary layer thickness, noise 
frequency, and duct dimensions. Since all of the correlated quantities men- 
tioned above involve only the optimum quantities, an off-optlm\jm estimate pro- 
cedure is also provided. This involves the approximate equation of reference 8 
in which the off-optimum behavior is shown to be uniquely determined by the op- 
timum impedance and damping along with the actual off-optimum wall impedance. 

The procedure outlined in this paper requires only the addition of the quantify- 
ing of the acoustic power-? distribution. 

Methods for estimating the acoustic power-? distribution from the far field 
directivity pattern are nearing completion and should be available soon. A more 
desirable method using direct duct measurements with flush-mounted wall pressure 
transducers is currently being studied. 

The problem of changes in acoustic power-? distribution in going from a 
hardwall duct to a soft wall section is discussed. 


SYMBOLS 


A function of eigenvalue phase 

angle, eq. (15) 
c speed of sound, m/sec 

D duct diameter, m 

AdB sound power attenuation, dB 

AdB maximum possible sound power 

attenuation, dB 

F boundary-layer refraction 

function, eq. (11) 

f frequency, Hz 

G function of maximum possible 

attenuation, eq. (18) 

J Bessel function of first kind, 

m , 

order m 

K wave number, kx 

k (i)/c, m“^ 


L acoustic liner length, m 

Mq axial steady flow Mach number, 

free-stream uniform value 
m spinning mode lobe number 

(circumferential order) 

N normalized expected number of 

modes versus cutoff ratio 
P acoustic pressure, N/m^ 

R amplitude of eigenvalue a 

Rjj^ hardwall eigenvalue 

r radial coordinate, m 

rQ circular duct radius, m 

t time, sec 

Vg group velocity (8w/3K) 

X axial coordinate, m , , 

X0 V 

a complex radial eigenvalue (ot = Re ) 
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6 AdB^/AdB 

6 boundary-layer thickness , m 

e dimensionles s boundary-layer 

thickness, S/rg 
optimum specifie acoustic 
impedance 

^mO optimum specific acoustic 

impedance for e = 0 
n frequency parameter, fD/c 

0 specific acoustic resistance 

0^ optimum specific acoustic 

resistance 


U radial mode number 

5 cutoff ratio 

5hW cutoff ratio in hardwall duct 
0 attenuation coefficient 

T propagation coefficient 

^ angular coordinate, rad 

(j> phase angle of eigenvalue, deg 

X specific acoustic reactance 

Xm optimum specific acoustic re- 

actance 

0 ) circular frequency, rad/sec 


DEFINITION OF THE CUTOFF RATIO 


Some preliminary expressions are given here to establish the notation and 
terminology. The modal pressure solutions are given by (in the uniform flow 
region) , 

P = J (1 

\ro/ 

where P, a, 0 , and t should actually have tn, y subscripts to associate 
them with the m, y mode. For soft walled ducts the radial eigenvalue is com~ 
plex and is given by 

a = Re^*^ (2) 

The damping and propagation coefficients are given by 


or 


0 + it = 


-imq + i^i- (1 - 


1 



( 3 ) 


0 + ix 



+ i 



(1 - m2) 



1 



(cos 2(fi + i sin 2(j)) 


where n is the frequency parameter given by 

ri = fD/c 


( 4 ) 


( 5 ) 
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For hardwall duct the definition of cutoff ratio is quite direct since the 
eigenvalue a is real and can be given by (in a manner similar to ref. 9), 



( 6 ) 


This definition causes the expression in the radical of equation (4) (with 
(j) = 0) to change sign at Cuy “ 1 and causes the pressure to be damped with 
distance for g < 1. 

For soft walled ducts the definition of cutoff ratio is not so simple. 

The modes possess propagating characteristics at all frequencies so there is no 
precise cutoff. A definition of cutoff ratio used in reference 7 was 


g = 



TTTI 


- Mq) cos 2(j) 


(7) 


which causes the real part of the radical in equation (4) to be zero at 5 = 1. 
The definition was quite arbitrary with the only advantage being that it re- 
duced to the hardwall definition when (j> ^ 0. A better definition might be 

TTTi Vcos ^ - /3 sin 

5 ( 8 ) 

rVi - Mq 

which is thought to be a new result. This was derived by insuring that the 
group velocity (Vg = 3u)/3K) be at a minimum when g - 1, which implies a mini- 
mum acoustic power propagation. Fortunately there is not much difference be- 
tween the 5 definitions for the small angles tp encountered at the optimum 
impedance. For the largest difference C from equation (8) is about 0.87 of 
that from equation (7). Thus the calculated results of reference 7 are used 
here without modification. 


MODE CUTOFF RATIO AS THE BASIC PROPAGATION PARAMETER 


In this section the cutoff ratio will be shown to be the basic parameter 
governing noise propagation in acoustically lined ducts. This will be done by 
showing that the optimum wall impedance for all of the modes can be accurately 
correlated by the cutoff ratio alone and that the maximum attenuation at this 
optimum can be adequately correlated by the cutoff ratio. All of the calcula- 
tions presented here were obtained using the calculation procedure of reference 
10. A uniform flow region was assumed in the duct interior with a boundary 
layer present near the wall. The classic uniform-flow sound propagation solu- 
tions were coupled to a Runge-Kutta integration solution through the boundary 
layer. The definition of modal optimum impedance is the same as in references 
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2 and 10 as well as that of reference 11 but with the additional considerations 
of Mach number, boundary-layer thickness, and higher order modes* 


Optimxjm Wall Impedance 


The discovery that the optimum liner wall impedance depends only on mode 
cutoff ratio was documented in reference 7* Figure 1 is repeated from refer- 
ence 7 for completeness. Figure 1 shows sample optimum impedance calculations 
plotted in the wall impedance plane. The conditions used are those for a 
General Electric TF-34 engine. Note that a common locus of optima is evident 
despite the wide range of modes used. Only the first radial of each of the 
higher lobe number modes deviates from this common line, and this deviation is 
quite small. (The first radial is the furthest point toward the left side for 
a given lobe number.) Two coincident modes are singled out and compared in the 
insert table. The only thing these two modes have in common is the cutoff 
ratio. Additional results were shown in reference 7 which illustrated the ex- 
cellent correlation of optimum impedance with cutoff ratio. 


Maximum Possible Sound Power Attenuation 


The maximum possible sound power attenuation is generally expressed as 
AdB/(L/D) and plotted against the frequency parameter (n ~ fD/c) as in figure 5 
of reference 2. This type of plot has been recast in terms of the cutoff ratio 
(fig. 2). Several radial modes (y = 1, 2, 5, and 10) for lobe numbers of 
m = 1, 7, and 20 are shown. In each case for a given m, the y - 1 curve is 
the lowest and the attenuation increases monotonically with increasing y. Ex- 
cept for the first two radial modes of the lowest lobe number (m = 1), the 
curves cluster together. If an average curve is used through the cluster of 
curves, most of the modes will be adequately represented with the maximum error 
deviation from the average being about a factor of two for the lowest order 
modes. In a multimodal liner design, this error in only a few of the lower 
order modes is not anticipated to be significant. In those cases where a few 
low-order modes are known to carry the bulk of the acoustic power, the method 
proposed in this paper should not be used since for these conditions the direct 
modal approach is both simpler and more accurate. 

With the exception noted previously, the maximum attenuation of the multi- 
tude of modes can be adequately represented by a single function of the cutoff 
ratio alone. The equations involved with the attenuation will be given in the 
next section where approximate correlations are discussed. 


APPROXIMATE CORRELATING EQUATIONS 


In this section approximate equations are developed for the optimum imped- 
ance and the maximum possible attenuation. The correlation between exactly 
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calculated optimum impedance and cutoff ratio is considered to be firmly estab- 
lished, but the optimum impedance correlation equation given here must be con- 
sidered preliminary. If more exact results are required at this time, the com- 
plete calculation procedure of reference 10 is suggested with a representative 
mode used at each cutoff ratio value under consideration. 


Optimum Impedance Correlation Equation 


The correlating equation was derived using the approximate equation of ref- 
erence 10 as a starting point, which in turn was derived from the thin boundary- 
layer approximation theory of reference 12. The starting point from reference 
10 is 


(1 + e)C. 


mO 


’m 


1 - iFC 


mO 


6 + ix 

m ^m 


(9) 


where is the optimum impedance with a boundary layer and where 

^mO " ®m0 

Is the optimum impedance with the boundary-layer thickness e = 
quantity F is given by 


(10) 
0. The 


■■(=?)(■•?) 

where the first term is the simplest form of the equation in reference 10 and 
the second is an empirical correction needed in the vicinity of unity cutoff 
ratio. This correction is needed since the results of reference 10 and thus 
presiimably reference 12 are not valid near cutoff. 


The quantities in equation (10) must now be cast in terms of the cutoff 
ratio 5 if equation (9) is to be a function of cutoff ratio as it is known to 
be from the exact calculations. Because of limited space, the derivations can 
not be included here. Equations (21) to (23), (30) and (31) of reference 8 a- 
long with equation (31) of reference 2 were used. The variables Mg and q- 
on which boundary-layer refraction effects strongly depend (ref. 10) were car- 
ried intact through the derivation. Certain liberties were taken with the other 
variables such as replacing nearly first or second powers of the mode numbers 
with first and second powers of the eigenvalue. The eigenvalue had to be re- 
covered in the equations in order to introduce the cutoff ratio from equation 
(8). The final equation, which must be considered as empirical, is 


r ~ S 

’mO ~ (1 + Mq) 




1 » - V f 


( 12 ) 
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Equation (12) is surprisingly accurate for the zero boundary-layer thickness 
optimum resistance, but large percent errors can occur in the reactance for very 
small values of reactance. 

With equations (11) and (12) used in equation (9) the optimum wall imped- 
ance with a boundary layer can be calculated as a function of cutoff ratio. 

These approximate calculations are compared with the exact calculations (from 
ref. 7) in figures 3 and 4. The approximations predict the gross behavior of 
the exact calculations and are probably accurate enough for most liner design 
studies. 


Maximum Attenuation Correlation 


An expression for maximum possible sound power attenuation can be derived 
by using the real part a of equation (4) and the cutoff ratio from equation 
(8) in the following expression (ref. 2) 


AdB 


m 


L/D 


-l7.4iTria 


which then yields 
AdB 


m 


-17. 4R 


L/D 


^ 2(1 - M^) L 

where A is given by 


(a) " ^(t) (f) " (a) 


1/2 


A = ^cos - /j" 


sin 


2i 


(13) 


(14) 


(15) 


Equation (14) was used to generate the curves in figure 2 with the eigenvalues 
(R, (j>) used for each mode and with Mq “ 0. An averaged equation can be ob- 
tained by using the values of R and (j) near the center of the cluster of 
curves such as the 20, 1 (m, y) mode (R = 26.662, = 5.46°) or the 7, 10 mode 

(R = 41.881, 4 = 3.53°). An approximate form of equation (14) can be derived 
for large 5 as 


AdB 

m 

L/D 


-8.7RA sin 2(ti .. -40 

■ ■ ■ ■ 

sVl - Mq ?Vl - Mq 


(16) 


where the average-curve values of R and <t> were used to arrive at the final 
expression. Th is compa res favorably with the expression in reference 7 except 

that the term Vl - is missing in reference 7. When expressed on a modal 
basis, the Mach number would not appear in equation (16) as discussed in refer- 
ence 2, but it is reinserted along with the cutoff ratio when equation (8) is 
used. Equation (16) is valid for the linear portion of the curves in figure 2 
but it will underpredict the maximum attenuation near g = 1. 
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OFF-OPTIMUM ATTENUATION 


As seen in figure 1 for any chosen impedance, at best only one value of 
cutoff ratio would be at an optimum. Since a distribution of modes and their 
associated cutoff ratios would be the usual case encountered, off-optimum damp- 
ing must be considered. Since the optimum impedance and maximum damping are 
now known for any cutoff ratio value, the approximate attenuation equation of 
reference 8 is ideally suited for use here. This equation is expressed as 

9^ - 2e [6ni + G(S - 1)] + 

6 ( 3 ^ 

8(B + 1) 


0 ^( 28 ^ - 8 + 1 ) 

- 3(3 + 1 ) 


+ 2G(8 - 1) 


== 0 


(17) 


where 


G = 8.7 L/D 

SdB^d + 

and 


(18) 


AdB 

_ m 

■" AdB 


(19) 


The of f-optiniuin attenuation to be solved for is AdB occuring at an impedance 
given by 0 and x with optimiam input values ^ with the 

usual design inputs of L, D, and When operating off-optimiim AdB cannot 

exceed AdB , and g is always greater than unity. Thus a possible procedure 
for solving equation (17) is to increment g upward from unity until the equa- 
tion is satisfied and then solve for AdB from equation (19). 


OUTLINE OF USE OF THE ACOUSTIC LINER DESIGN PROCEDURE 


In the example that follows the use of the equations presented in the pre- 
ceding sections will be illustrated. The final element of the technique that 
is needed is the distribution of acoustic power as a function of cutoff ratio. 


Because of space limitation as well as 
of this subject, the equations will be 
for illustrative purposes only. 

The modal population density as a 

dC 

If equation (20) is integrated between 
modes between these limits is 


the preliminary state of the development 
presented without proof and are intended 

function of cutoff ratio is expressed as 

= (20) 
5 

and ^ 2 » normalized number of 
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(21) 



Note that N is nonnalized since, if ^ “ "» then N = 1. For 

simplicity equal energy per mode will be used, which may be a reasonable as- 
sumption for static test data (ref. 5). Then the acoustic power between 5-i 
and also given by equation (21). Next, choose 10 increments so that 0.1 

of the power falls in each interval with 0.05 of the power on each side of the 
interval center. The 10 g intervals thus have centers located at 5 - 1.026, 
1.085, 1.155, 1.24, 1.35, 1.49, 1.69, 2, 2.58, and 4.47. 

Now the attenuation calculation can be made for each value of cutoff ratio 
5 at the desired values of resistance 9 and reactance x* Assume that 
boundary-layer thickness e, Mach number Mq, frequency parameter n, and the 
allowable duct diameter D and length L are known from other considerations. 
For each center value of g, then, calculate the optimum impedance components 
0JU and Xm from equation (9) using equations (11) and (12). Next, calculate 
the maximum possible attenuation AdBjjj for each C from equation (14) using 
equation (15) and the R and 4 values given just after equation (15). Now 
all the inputs are available (6, x» 0m> Xm» calculate 6 from equa- 

tion (17) and then AdB from equation (19), again, for each of the 10 values of 
5. The estimated liner output acoustic power can then be calculated, and an 
overall AdB calculated by applying the AdB for each g catagory to its re- 
spective input power, summing the output powers, and comparing this sum with the 
total input power. The calculation is now complete at the selected value of .6 
and X* If ^ multimodal optimization study is being made, new values of 9 and 
X would be selected, and the calculations repeated until the total attenuation 
is maximized. 


CONCLUDING REMARKS 


The acoustic liner evaluation method presented in this paper should provide 
a useful alternative to the more usual modal analysis approach. Some of the 
problems in the modal approach, which are not problems in the present approach, 
are as follows: The phase speed of a mode is inversely proportional to the pro- 

pagation coefficient t given by equation ( 3). If t he cutoff ratio from equa- 
tion (6) is used in equation (3) , the term A - 1 /%^ will be found to contain 
all of the modal information. For well propagating modes (g >> 1) this radical 
is essentially equal to one and the axial phase speed of all these modes is 
nearly identical. Also in a multimodal situation several modes may have almost 
equal cutoff ratios even though g w 1. These modes would also be indistin- 
guishable in an axial direction since they have the same axial phase velocity. 
Thus the modal acoustic power can not be uniquely determined by using axial 
microphone traverses, and radial traverses (which are undesirable) must be used. 
Modes with nearly coincident cutoff ratios do not present a problem to the 
method of this paper since they all behave similarly in the acoustic liner (with 
respect to optimum impedance and maximum attenuation) and are thus lumped to- 
gether . 
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Some simplifying conditions have been tacitly assumed to hold in the devel- 
opment of the technique presented here. It has been asstaned that modal cross 
coupling is not important in the attenuation calculations in the lined duct sec- 
tion; that is, the reduction of acoustic power of each mode can be calculated 
independently of all the other modes. Theoretically, the cross coupling should 
be considered (refs. 13 and 14), but for practical purposes this coupling might 
be negligible (0.5-dB error in the results of ref. 15). Although not directly 
affecting the results presented here, the problem of a possible change in acous- 
tic power-? distribution in going from a hard walled section to a soft walled 
section must be recognized. Could a power density distribution determined in a 
hard duct be used for attenuation calculations in the lined duct section or what 
modifications must be made? This can not be answered definitely at this time, 
but some insight can be offered. The cutoff ratio should be related to the 
angle of incidence with the duct wall. If modes are available in the soft duct 
with angles similar to those in the hard duct, then these modes should be ex- 
cited. Thus angle of incidence is preserved in much the same way as the lobe 
numbers are preserved. This analogy is not exact since modes with exactly the 
same lobe number are available in both duct sections while angles of incidence 
can only be approximately the same. Thus some scattering of acoustic power 
among the various cutoff ratios should occur, but in a multimodal situation this 
is not suspected to be an extremely important effect. A notable singular ex- 
ception occurs when a plane wave in the hard duct reaches the soft walled sec- 
tion. The plane wave with ? = » is scattered, into several soft wall modes 
with finite and possibly even small cutoff ratios (depending on frequency param- 
eter). This situation does not comply with the conditions assumed earlier, 
since there is no mode in a very soft duct that matches the angle of incidence 
of the plane wave. Perhaps the most serious assumption of all is that an 
acoustic power-? distribution is more available than an acoustic power-modal 
distribution. In reference 7 it was implied that the far-field directivity pat- 
tern is intimately related to the acoustic power-? distribution. This approach, 
which has been pursued and is nearing completion, should allow at least a crude 
approximation to the power distribution. Also, the present method offers the 
potential for avoiding some of the problems associated with modal measurement. 
Measurements of the acoustic power-? distribution using wall mounted micro- 
phones in the duct should be developable and ultimately available. 
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Figure 1. - Example higher order spinning mode optimum im- 
pedance Ixus. Frequency, 2890 hertz; frequency parameter, 
9. 47; Mach number, -a 36; boundary layer, 6/rQ = a 059. 


figure 2. - Maximum possible attenuation as 
function of mode cutoff ratio. Zero Mach 
number. 
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Figure 3. - Resi stance correlation compared with exact cal- 
culations. 


Figure 4. - Reactance correlation compared with exact cal- 
culations. 
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